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ABSTRACT 

Using the theoretical model of S. Yoksan (Solid State Commun. 
78, 436 (1991),we have studied the isotope effect in high T c 
superconductors. This model takes into account the nearest 
neighbor interlayer hopping. Our theoretical results indicate that 
the value of isotope exponent a of high T c superconductor is very 
small. 
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INTRODUCTION 

One of the most interesting 
experimental properties of high T c 
superconductors is the small isotope effect. 
The suppression of this effect is not 
explained with the conventional BCS theory 
predicts that the critical temperature T c and 
the isotope mass M are related by T c M~ a , 
where a = 0.5 for all elements. 1 " 3 

To explain the almost absence of the 
isotope effect in high T c compounds, many 
authors have proposed various mechanisms 
including resonating valence bond, 4 
excitons, 5 plasmons 6 and ant ferromagnetic 
spin fluctuation mediated pairing. 7 Fund and 



Kwok suggested a local volume change due 
to oxygen deficiency as a cause of this 
effect. Su et.al. interpreted experimental 
data as evidence for the free carrier negative 
U-centre interaction, Seki and Tanabe 10 
investigated the influence of the singularity 
enhanced density of states in a bond 
asymmetric high T c material. Daemen and 
Overhauser 11 surprisingly found that the 
existence of a short range attraction in 
addition to the conventional phonon pairing 
interaction suppresses the isotope effect 
significantly at high temperatures. 

In these calculations the authors 
assumed the constant electronic density of 
states at the Fermi level. As is well known, 



Journal of Pure Applied and Industrial Physics Vol.2, Issue 2, 1 April, 2012, Pages (84-141) 



99 Bharat Singh, et al, J. Pure Appl. & Ind. Phys. Vol.2 (2), 98-103 (2012) 



high temperature superconductors have a 
two dimensional structure. Various 
calculations 12 17 shows that there is a 
logarithmic singularity in density of states, 
one therefore expects that this singularity 
plays a dominant role in causing the 
anomalous behavior of the isotope effect. 
Theoretical Formulation 

Using the theoretical model given 
by Daemen and Overhauser 11 one assumes 
that the pairing interaction written as 



V kk ~ \_ v _ 



~Vi -V 2 JorO<E ha) D 



7 2 for hd) D < E < 4t ^ 

is independent of the phonon frequency. 
Here V! is the phonon mediated interaction 
and V 2 is the extra interaction which may be 
of non electron phonon origin. E is the 
electron energy measured from the Fermi 
energy, © D is the Debye cut-off energy and 
4t is the energy cut-off for V 2 . 

As it has been shown by experiments 
that in the oxide superconductors the carriers 
is still the Copper pair, one thereby assumes 
that a solution of the BCS integral equation 



A(fe) 



_ 1 ^ VQc - g)A (q) ^ V f 2 + A (q) 



2 V V^ 2 + A 2 (?) 



tanh- 



2k p T 



(2) 

Based on equation (1) has the following 
form 



_ (A 0 (r)A 1; if E < ha) D 
(T)A 2 , if E > hco D 



(3) 



Here A-l and A 2 are temperature independent 
constants. Upon substituting equation (3) 
into (2), one obtains the following equations 



1- Ci + V 2 ) JV) 



= A 2 V 2 E 2 (T) 



A^E^T) = A 2 [l- V 2 E 2 (T)] 



(4) 



when 



Elm - s; 



hto D dEN (E) 



tanh 



E 2 + Aq &l 



dEN (E) 



tanh 



E 2 +A 2 a A 2 2 



E 2 +Aq a\ 



2k p T 



(5) 



One obtains equation (5) by replacing the 
summation over q by the energy integrals. 
Now in the Hubbard model on a two- 
dimensional square lattice, the Hamiltonian 



is 



//= ^tC+ + C io + 



(6) 



IjO 



here t denotes the transfer integral, and U the 
renormalized on site Coulomb interaction, 
the number operation n io — C^C io . In case 
of the nearest neighbor hopping, the band 
energy is E (k)= -2t (cosk x + cosk y ).This 
gives rise to the density of states in the form 



N(E) = 



2* 2t 
0 



| E | < 4t 
|E| > 4t 



(7) 



When K is the complete elliptic integral of 
the first kind, and N (E) is non-zero only 
when 0 < E < 4t, for E 0, N (E) takes the 
form 



N (E) = — — In 

v ' 2n 2 t 



I E I 

Il6tl 



(8) 



AT T = T c , A 0 (T) = 0, with the aid of the 
expression for N(E) one can be see that the 
integrals in equation (5) are dominate around 
E=0, so one approximate equation (5) as 

Z (7;) = 2nH F (2^) 
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With the function F defined by 
F (zfc^) = S 6 



hco ° dx In x / 8tx 
tallh [ Wc 



(10) 



When equation (4) are compatible one 
arrives at the following condition 



2k p T c 



2k p T c 



2k„T r 



At 



- F 



hai r , 



= 0 
(11) 



Here one has introduced the variable 
X = VJ2Ti 2 t and o = V 2 /27i 2 t . 



T c CALCULATION 



To calculate 13 the critical temper atureT c , one 
evaluates F and finds that 



(k js T c 



\2kpT c ) 2 \ 4t ) \ 8t 

Putting F in equation (11) and rearranging, 
one arrives at the equation for T c . 



0.6646 



k p T c = 1.13ha) D expl^t-^ — — 

p I ln(ftw D /16t) 

Where 



(13) 



(14) 



^ 2 ^ n (ft(D D ) ^ n (ftdl D ) 

In the pure electron-phonon mechanism, i.e. 



when V-l 0 and V 2 = 0 and o = 0 
equation (13) gives. 



fc^ = 1.13ft&j D expj j^ Vw 

Uir 2 t " V I6t J J 



(15a) 



In the dominate non electron-phonon 
process, i.e. when V 2 > > V 1; we have 
X -> 0 and o 0 equation (13) gives. 



fc^r c = 8texp - 



47T 2 t 



(15*) 



Here the quantity of l/2n 2 t\ln(hOi) D / 
16t) \ is analogous to the density of states at 
the Fermi level of the BCS theory. If this 
quantity is large, r c ~<y D .Thus the 
logarithmic singularity yields a function 
|l/ln(/iCL) D /16t)| giving an important 
enhancement of T c . 

For V 2 0, one calculates T, as a 
function of a/1 for different ratios hu> D /16t. 

a CALCULATION 

The isotope effect exponent can be derived 
from equation (13) as 



a = 



d lnT c 
d In M 



(16) 



In the phonon-mediated pairing 
scheme cd D is proportional to M~ 1/2 , when 
M is the average mass of the atom 
constituting the specimen under 
consideration, using equation (13) in 
equation (16), one finds 



[1 - X (0.6646 + In ^ In 

V 16t 1.134fta) D y 



V /1.134ft(D D/ ) 
ln(ftoj D /16t) 



(17) 
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in the limit of low T c , it is straight forward to show that equation (17) can be 
estimated as 



a= ifl-h-^lnf^xln^^-)! 2 ) 

2 ( L V 16t 1.134fi(<) D /J J 

Here the quantity 



-4- In 

2n 2 t 



(18) 



(19) 



in this model is equivalent to NoV^ in equation (12). 

In a refined version of a, one calculates it directly by differentianting equation (11) with 
respect to M and get 



H(A(u D /fe^r c ) 



[1 + X{ha> D /k p T c )] 2 



P'cJ 

ln 4 tanh 2fe^ 



H (z&)] 



2 +X.F 
-l 



»fe)'-(^fc)- H (^fc)' 

ry/ 16t dx /8tx\ 

— tanh — - 

J 0 X V fc /? r c/ 



Table T, 



High T c Superconductor 


T C (K) 


®D (K) 


t(K) 




a theo 


a expt 


La 2 _ x Sr x Cu0 4 


37 


240 


5000 


0.02 


0.19 


0.16 


YBa 2 Cu 3 0 7 


92 


450 


6000 


0.01 


0.09 


0.018 
±0.004 


Bi 2 Sr 2 Ca 2 Cu 3 O 10 


107 


535 


6925 


0.025 


0.048 


0.026 
±0.02 


Bao.6Ko.4BO3 


32 


186 


4783 


0.022 


0.40 


0.042 
±0.05 


(Tl 0 . 5 Pb 0 . 5 )Sr 2 Cu 2 Cu 3 O 9 


122 


639 


7447 


0.039 


0.037 


0.027 
±0.003 



DISCUSSION OF RESULTS within the conventional BCS theory. By 

passing that there exist the usual attractive 
In this paper, we have studied the interaction, Originating from the phonon 
isotope effect in high T c superconductor exchange and an attractive interaction from 
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the unknown origin. One calculates the 
transition temperature T c and the isotope 
effect exponent a. The possible singular 
behaviour in the density of states and the 
relative importance of the attractive 
interaction T c and a are investigated. It is 
found that the singularity contributes to the 
enhancement of T c . For high critical 
temperature, the second attractive 
interactions must be the stronger of the two. 
One also observes that the behaviour of the 
exponent a depends sensitively on the 
relative magnitude of the two interactions. 
The results are shown in table 7\. 
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